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O Explosive nuclear burning in astrophysical environments produces unstable nuclei 

3 , which again can be targets for subsequent reactions. In addition, it involves a 

£Z< ' large number of stable nuclides which are not fully explored by experiments, yet. 

Thus, it is necessary to be able to predict reaction cross sections and thermonu- 
clear rates with the aid of theoretical models. Such predictions are also of interest 
, for investigations at radioactive ion beam facilities. An extended library of theo- 

> I ■ retical cross sections and reaction rates is presented. The problem of a+nucleus 

potentials is addressed and new parametrizations presented. The problem of prop- 
erly predicting cross sections at low level densities is illustrated by the 62 Ni(n,7) 
reaction. 



1. Introduction 

The majority of reactions in astrophysics involving the strong interaction 
can be described in the statistical model 1 . In predictions of cross sections 
for astrophysical applications, slightly different points are emphasized than 
in pure nuclear physics investigations. Firstly, one is confined to the very 
low energy region, from thermal energies up to a few MeV. Secondly, since 
most of the ingredients for the calculations are experimentally undeter- 
mined, one has to develop reliable phenomenological or microscopic models 
to predict these properties with an acceptable accuracy across the nuclear 
chart. Therein one has to be satisfied with a more limited accuracy as 

1 
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compared to usual nuclear physics standards. Considering the substan- 
tially larger uncertainties in many astrophysical scenarios, this seems to be 
adequate. 



2. Statistical model calculations 

A recently published large-scale reaction rate library includes neutron-, 
proton-, and a-induced reactions on all target nuclei from Ne up to Bi from 
proton-dripline to neutron-dripline 2 ' 3 . Due to the fact that many very 
short-lived nuclides can be produced in astrophysical sites, it is necessary 
to provide cross sections and rates for about 4600 targets and 32000 reac- 
tions. These numbers show that theory will always play a major role in 
providing cross section, despite the potential of future Rare Isotope Accel- 
erators. The calculations were performed with the Hauser-Feshbach code 
NON-SMOKER 4 which is especially tuned to such large-scale predictions. 
Details of the nuclear properties used are given elsewhere 2 . This rate set 
has already be adopted as a standard for nucleosynthesis in stellar evolution 
and in type II supernovae 5 . 

Fits to the astrophysical reaction rates - ready for direct astrophysical 
application - as well as tables of cross sections, reaction rates, and nuclear 
inputs for all possible reactions with light projectiles can also be downloaded 
from http://nucastro.org/reaclib.html . 



3. Optical a-nucleus potentials 

There have only been few attempts to derive global optical potentials for 
a-projectiles 6 and most of them are only valid at a-energies larger than 
30 McV. Due to the high Coulomb barrier and nuclear structure effects 
defining the imaginary part of the potential it is difficult to obtain a global 
potential at astrophysical energies. Elastic a-scattering data can constrain 
the real part of the potential 7,8 and detailed analysis can also improve 
on the imaginary part 9 ' 10 , describing the absorption into other channels 
than the elastic scattering, i.e. the Hauser-Feshbach channel. Due to the 
scarcity of data for intermediate and heavy nuclei, attempts to improve on 
the potential are mostly concentrating on single reactions 14 ' 10 . More global 
approaches suffer from the lack of data to confine their parameters 6,11 . 

We have tried to find a potential for the A = 140 mass region by si- 
multaneously fitting data for 143 Nd(n,a) 140 Ce 12 , 147 Sm(n,a) 144 Nd 13 , and 
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Figure 1. Cross sections, reaction rates, and S-factors from a simultaneous x 2 fit of the 
Fermi-type energy-dependent a+nucleus optical potentials of three reactions (see text). 
The dashed lines are the statistical model calculation. The errors on the 147 Sm(n,«) 
rates were assumed to be 10%. 



144 Sm(a,7) 148 Gd 14 . The optical potential is parametrized as 

V . W(E) 



V(r,E) 



l + exp(^^) l + exp(r^i) 



(1) 



Different parameters for the potential geometry and the energy depen- 
dence of the depth of the imaginary part were explored 15 . We did not 
find significant differences between using a Brown-Rho shape 8 W{E) = 
W {{E - E a ) 2 )/((E - E ) 2 + A 2 ) or a Fermi-type shape 14 W(E) = 
Wo/(l + exp((i?* - E)/a*)) of the energy dependence. For the latter we 
found the parameters E* = 18.74 MeV, a* = 2.1 MeV, with all other pa- 
rameters as in the previous paper 14 . The Brown-Rho best fit was obtained 
with E — 6.35 MeV and A = 28.4 MeV, with the same fixed parameters 
V = 162 MeV, r r = 1.27 fm, a r = 0.48 fm, W a = 19 MeV, r v = 1.57 fm, 
ay = 0.6 fm. The results from the simultaneous fit of three reactions are 
shown in Figs. 1 and 2. 

Despite the fact that the considered targets are in the same mass re- 
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Figure 2. Cross sections, reaction rates, and S-factors from a simultaneous x 2 fit of the 
Brown- Rho energy-dependent o+nuclcus optical potentials of three reactions (see text). 
The dashed lines are the statistical model calculation. The errors on the 147 Sm(n,o) 
rates were assumed to be 10%. 



gion, the derived parameters also describe acceptably well the reaction 
96 Ru(a,7) 100 Pd 16 . However, it is remarkable that even better overall 
agreement with all four reactions can be obtained when using a mass- and 
energy-independent potential of Saxon- Woods form for the real and imagi- 
nary parts (see Fig. 3). The real parameters are given by V = 162.3 MeV, 
r r = 1.27 fm, a r = 0.48 fm, the imaginary ones by Wq(E) — Wv = 25 
MeV, rv = 1-4 fm, ay = 0.52 fm. Thus, the real part is identical to the 
potential by Somorjai et al. but without energy dependence, whereas the 
imaginary part is similar to the one used in McFadden & Satchler 14 . Since 
the McFadden & Satchler parameters were derived from extensive elastic 
scattering data it seems reasonable that they are applicable to a wider range 
of targets. The Somorjai et al. parameters were derived for one reaction 
only but seem to work also for the nuclides investigated here. Certainly, at 
very low a-energies an additional energy-dependence has to be introduced. 
Here, we do not show our results from fitting each reaction separately. Ob- 
viously, potentials fitted to single reactions can describe those - but only 
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Figure 3. Results for four different reactions using the energy-independent potential (sec 
text). The dashed lines are the statistical model calculation. Note that the 147 Sm(n,a) 
result is rcnormalized by a factor 1/1.4. 



those - even better. 



4. The 62 Ni(n,7) case 

For neutron-induced reactions at low energies, close to magic numbers, 
and far off stability where low separation energies are encountered, another 
problem emerges. In such targets, the level density is too low to allow 
the application of the statistical model 1 . Also for other nuclides it is 
not straightforward to bridge the region of thermal energies to the region 
of overlapping resonances where the Hauser-Feshbach formalism can be 
used. Single resonances and direct reactions become important. This is 
also an issue for neutron-rich nuclei in the r-process path with low neutron- 
separation energies. 

As an example for the difficulties in extrapolating thermal data to s- 
process energies of up to a few hundred keV, the reaction 62 Ni(n,7) 63 Ni 
is discussed here. Two compilations give disagreeing 30 keV cross sec- 
tions 17 ' 18 J based on the same thermal data. Both extrapolations assume 
s-wave behavior of a direct capture component. The more recent one in- 
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Figure 4. Direct neutron capture Maxwellian averaged cross section of 62 Ni. The final 
value is given by adding the resonant contribution to the "total" direct term. Upper and 
lower limits on the direct components are from experimental errors on the input, i. e. in 
the thermal scattering length and the spectroscopic factors. 

eludes a sub-threshold resonance contributing to the thermal cross section. 

We have calculated the direct capture component using DWBA and 
found a considerable p-wave contribution which enhances the cross section 
at 30 keV 19 . Thus, even when including the subthreshold resonance, the 
30 keV value is coincidcntally similar to the value in the older compilation 
(Fig. 4). However, also the general energy dependence of the cross section 
is altered. Resonances were also included but they only contribute less than 
15%. The enhanced cross section has an important impact on s-processing 
in massive stars. A previously seen overproduction of 62 Ni in stellar models 
can be cured when using our enhanced rate because of increased destruction 
of this nucleus with the larger neutron capture rate 5 ' 19 . 

5. Conclusions 

Despite considerable successes in the prediction of cross sections and reac- 
tion rates close to and far off stability, the description of certain nuclear 
inputs, such as optical a-potentials, still needs to be improved. It is also 
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still unclear whether nuclear properties far off stability can be predicted 
with sufficiently high accuracy. Although future advances in microscopic 
theories may alleviate that problem, experimental data is clearly needed. 
Rare Isotope Accelerators will make it possible to study highly unstable 
nuclides but also "classical" nuclear physics experiments with stable or 
long-lived nuclei are indispensable. They can provide the systematics for 
global descriptions and shed light on the interaction of different reaction 
mechanisms. 
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